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Abstract 

A computation method of algebraic local cohomology with parameters, associated with zero¬ 
dimensional ideal with parameter, is introduced. This computation method gives us in particular 
a decomposition of the parameter space depending on the structure of algebraic local cohomology 
classes. This decomposition informs us several properties of input ideals and the output of our 
algorithm completely describes the multiplicity structure of input ideals. An efficient algorithm 
for computing a parametric standard basis of a given zero-dimensional ideal, with respect to an 
arbitrary local term order, is also described as an application of the computation method. The 
algorithm can always output “reduced” standard basis of a given zero-dimensional ideal, even 
if the zero-dimensional ideal has parameters. 
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1. Introduction 


Local cohomology was introduced by A. Grothendieck in (|Grothendiecld . 119671) . Sub¬ 
sequent development to a great extent has been motivated by Grothendieck’s ideas 
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Brodmann. M. P. and Sharp. R. ^ ( 1998f l: Lvubeznik. G. ( 2002h . Nowadays, local co¬ 
homology is a key ingredient in algebraic geometry, commutative algebra, topology and 
D-mo dules, and is a funda mental tool for applications in several fields. 

In ( Taiima et al.l . [2009ll . we proposed, with Y. Nakamura, an algorithmic method to 
compute algebraic local cohomology classes, supported at a point, associated with a given 
zero-dimensional ideal. We described therein an efficient method for computing standard 
bases of zero-dimensional ideals, that utilize algebraic local cohomology classes. The un¬ 
derlying idea of the proposed method comes from the fact that algebraic local cohomol¬ 
ogy classes can completely describe the multiplicity structure of a zero-dimensional ideal 
via the Grothendieck local du ality theorem. More recently in our result of ISSAC2014 
(INabeshima and Taiimal . 1201411 . we considered the Jacobi ideal, with deformation param¬ 
eter, of a semi-quasiho mogeneous hypersur face isolated singularity. By adopting the same 
approach presented in iTaiima et al. ( 20091) , we constructed an algorithm for computing 
algebraic local cohomology classes, with parameters, that are annihilated by the Jacobi 
ideal. As an application, we obtained a new method to compute parametric standard 
bases of Jacobi ideals associated with a deformation of semi-quasihomogeneous hyper¬ 
surface isolated singularities. 

In this paper, we address the problem of finding an effective method to treat algebraic 
local cohomology classes with parameters associated with a given zero-dimensional ideal 
with parameters, that works in general cases. 

In order to state precisely the problem, let X be an open neighborhood of the origin O 
of the n-dimensional complex space C" with coordinates x = {xi,X 2 ^. ■., Xn)- We assume 
that a set F of p polynomials fi, f 2 , ■ ■ ■, fp in (C[ti,..., <m])[a^] satisfying generically 
{a € A|/i(a) = ••• = fp{a) = 0} = {O} are given where are parameters. 

Let Hp he a set of algebraic local cohomology classes supported at the origin that are 
annihilated by the ideal generated by F. Then Hp is a finite-dimensional vector space if 
and only if the ideal (F) generated by F is zero-dimensional in the rings of formal power 
series. In such cases, there is a possibility that {a e A|/i(a) = • • • = fp{a) = 0} ^ {O} 
(the same meaning is that {F) is not zero-dimensional) for some values of parameters, 
because of parameters. As our aim is to construct algorithms for studying the structure of 
Hp and the multiplicity structure of (F) on A, it is necessary, beforehand if possible, to 
detect these values of parameters, that constitute constructible sets, from the parameter 
space for computing algebraic local cohomology classes. 

In the first part of this paper, we introduce a new notion of parametric local cohomol¬ 
ogy system as an analogue of comprehensive system to deal with parametric problems. 
We describe a new effective method to compute parametric local cohomology systems. 
The resulting algorithms compute in particular a suitable decomposition of parameter 
space to a finite set of constructible sets according to the structure of algebraic local 
cohomology classes in question. The key of the algorithm for decomposing is the use of 
a comprehensive Grobner bases computation in a polynomial ring with parameters. The 
algorithms for computing bases of Hp, is designed as dynamic algorithm in consideration 
of computational efficiency. The output of our algorithm, has the abundant information 
of the input ideal and provides a complete description of the multiplicity structures of 
parametric zero-dimensional ideals. 

In the second part of this paper, we describe algorithms for computing parametric 
standard bases as an application of parametric local cohomology systems. We show that 
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the use of algebraic local cohomology provides an efficient algorithm for computing stan¬ 
dard bases. Furthermore, the use of algebraic local cohomology transforms a standard 
basis of a dimensional ideal (F) with respect to any given local term order into a standard 
basis with respect to any other ordering, without computing the standard basis, again. 
In general, the computation complexity of standard bases, is strongly influenced by the 
term order, like Grdbner bases computation. Thus, this property is useful to compute a 
standard basis. 

Especially, our algorithm can output always “reduced” standard basis of a given zero¬ 
dimensional ideal, even if F has para meters. Note that, an alg orithm implemented in 
the computer algebra system Singular ( Decker. W. et all 2012 ) that compute standard 
bases does not enjoy this property. Moreover, in general, comprehensive Grobner ba- 


19921) in a polynomial ring does not have this 


sis (iNabeshimal . I2012t iWeispfenning. V, 
property, too. 

As we mentioned above, there are several applications of algebraic local cohomology. 
For examples, our algorit hm can be used to analyze properties of singula rities and defor¬ 
mations of Artin algebra ( larrobino and EmsalemLll978 : Iarrobin3.ll984l) . It is a powerful 
tool to study several problems relevant to zero-dimensional ideals. 

All algorithms in this paper, have been implemented in the computer algebra system 
Risa/Asir ( Noro and Takeshimal . Il992 ). 

This paper is organized as follows. Section 2 briefly reviews algebraic local cohomology, 
and gives notations and definitions used in this paper. Section 3 is the discussion of the 
new algorithm for algebraic local cohomology classes with parameters. This section is the 
main part of this paper. Section 4 gives algorithms for computing parametric standard 
bases for a given zero-dimensional ideals. 


2. Preliminaries 

In this section, first we briefly review algebraic local cohomology. Second, we intro¬ 
duce a term order for computing algebraic local cohomology classes and algebraically 
constructible sets, which will be exploited several times in this paper. Throughout this 
paper, we use the notation x as the abbreviation of n variables cci,... ,x„. The set of 
natural number N includes zero. K is the field of rational numbers Q or the field of 
complex numbers C. 

2.1. Algebraic local cohomology 

Let F[^q^{K[x\) denote the set of algebraic local cohomology classes supported at the 
origin O with coefficients in AT, defined by 

H[o]{K[x\) := yu^Ext'^^^^{K[x\/{xi,X 2 , .■,Xn)'",K[x\) 

where {xi,X 2 , ■.. ,Xn) is the maximal ideal generated by xi,X 2 , •. • ,Xn. 

Let A be a neighborhood of the origin O of AT". Gonsider the pair (X,X — O) and 
its relative Cech covering. Then, any section of (Aja;]) can be represented as an 
element of relative Cech cohomology. We use the notation ^ c\ for representing 
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an algebraic local cohomology class in H^q^{K[x\) where c\ S K, 

with A = (Ai, A 2 ,..., A„) G N". Note that the multiplication is defined as 


r A+1 


A-j-1 — oc 


■■= < 


Xi ^ O-iy I - 1, . . . , Tly 


otherwise, 


where a = (ai,..., an) € N” and A + 1 — a = (Ai + 1 — ai, . .., A„ + 1 — an)- 

We represent an algebraic local cohomology class [y^] ^ polynomial in n 

variables X caX to manipulate algebraic local cohomology classes efficiently (on com¬ 
puter), where ^ is the abbreviation of n variables 5i, ^ 2 , ■ • ■, We call this representation 


"polynomial representation”. For example, let 7 /’ = 


be an algebraic 


local cohomology class where xi,X 2 are variables. Then, the polynomial representation 
of Ip, is 4^1 Ci + where variables ^ 1,^2 are corresponding to variables xi,X 2 - That 

is, we have the following table for n variables: 


Cech representation polynomial representation 


1 


Ai-l-l A 2 -I-I .^-^n+1 

U/1 2 ^Tl 


where ca G K. The multiplication for polynomial representation is defined as follows: 




a:“ * £^ := 


Xi ^ ai, % — 1 , . . . , Tly 


0, otherwise. 


where a = (oi,..., «„) G N", A = (Ai,..., A„) G N", and A —a = (Ai —oi,..., A„ — an)- 
We use “ * ” for polynomial representation. 

After here, we adapt polynomial representation to represent an algebraic local coho¬ 
mology class. We use mainly the following term order to compute algebraic local coho¬ 
mology classes. 


Definition 1. For two multi-indices A = (Ai, A 2 ,..., A„) and A' = (A(, A 2 ,..., A(j) in 
N”, we denote -< or A' ^ A if |£^'| < |£^|, or if = |£''| and there exists i,j G N 
so that A' = Xi for i < j and A' < Xj where |£^| = Xr=i general, this term order 

is called a total degree lexicographic term order. 

For a given algebraic local cohomology class ip of the form, ip = ca£^ + Xa'^a ca'£^ , 
Ca 7 ^ 0, we call the head term and , A' ^ A the lower terms. We denote the head 
term of a cohomology class ip by htpip). 

2-2- Strata and specialization 

We use the notation t as the abbreviation of m variables ti,... ,tm- (One can also 
regard t as parameters.) Let K be an algebraic closure field of K. For gi,--.,gq G K[t], 
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9q, i-e., Y{gi,...,gq) := {a G 


V( 5 i, ... ,gq) C K'^ denotes the affine variety of t/i, 

K"^\ 5 i(a) = ... = gg{a) = 0} and V(0) := 

We use an algebraically constructible set that has a form Y{gi, 

K™- where gi, ■ ■ ■ ,gq, g'l, - ■ ■ ,g'qi G K\t]. We call the form ¥( 51 , 
stratum. (Notation A, A', Ai,..., A; are frequently used to represent strata.) 

When we treat with systems of parametric equations, the n it is necessary to check con- 


,9q)\ a 


sistency of their parametric consiste nts. In several papers (|KaDur et al.L 1201(1 iMontes 


I 2 OO 2 I : [Suzuki. A. and Sato. Y.I . 2003 ). algorithms for checking consistency have been al¬ 
ready introduced. Thus, it is possible to decide whether Y{Qi)\Y{Q 2 ) is an empty set 
or not, by these algorithms where (5i, Q 2 C A'[t]. The details are in the papers. 

We define the localization of K\t] w.r.t. a stratum A C K'^ as follows: K[t]A = {f I 
c,b £ K[t],b{t) 7 ^ 0 for t G A}. Then, for every a G A, we can define the canonical 
specialization homomorphism as. '■ A'[t]A[a:] —>■ K[x\ (or as '■ A'[t]A[C] ^ When we 

say that as{h) makes sense for h G K{t)[x\, it has to be understood that h G if[t]A[a:] 
for some A with a G A^. We can regard as as substituting a into m variables t. 


3. Algebraic local cohomology with Parameters 

Let us assume that a set F of p polynomials /i, / 2 ,..., /p in (Ar[t])[a;] satisfying gener- 
ically {o G X\fi{a) = ■ ■ ■ = fp{a) = 0} = {O} are given where A is a neighborhood of 
the origin O of A”. Here, we regard t as parameters, and x, ^ are the main variables. 

We define a set Ap = ^seK’^^^asiF) fo be the set of algebraic local cohomology classes 
in A[^] that are annihilated by the ideal generated by F, where 

= {^ G K[^] I cra(/l) *i’ = 0-a(/2) * ijj = ■ ■ ■ = as{fp) * Ip = Q}. 

The ideal {F) at a G A"* is a zero-dimensional ideal if and only if is a finite¬ 

dimensional vector space. In this section we describe an algorithm for computing bases 
of the vector space Hp. More precisely, we describe algorithms for computing parametric 
local cohomology systems (see Definition 5 in this section). 

The new algorithm consists of the following three parts. 

(1) Decompose the parameter space A™ into safe strata and danger strata. 

(2) Compute bases of the vector space Hp on safe strata. 

(3) Compute bases of the vector space Hp on danger strata. 

3.1. An algorithm for testing dimensions of a parametric ideal 

Since polynomials /i,/ 2 ,...,/p have parameters, there is a possibility that {a G 
X\fi{a) = ■ ■ ■ = fp{a) = 0} 7 ^ {O}. As our aim is to construct algorithms for studying 
the system F on A, it is necessary, beforehand, to take away these values of param¬ 
eters that constitute constructible sets from the parameter space for computing local 
cohomology. 

Here, we describe an algorithm for decomposing A™ into S = {Ai,...,Afc} and 
T) = {Afe_|_i,..., A/} where {F) is zero-dimensional on A^ and nonzero-dimensional on 
Aj in a polynomial ring, for 1 < i < k,k + 1 < j < 1. This decomposition is possible 
by mainly computing a comprehensive Grobner system of F. We adopt the following 
definition of comprehensive Grobner systems, because this definition is suitable to com¬ 
pute dimensions of ideals in the algorithm ZeroDimension. (The following definition is 
different from the original one). 
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For any g € and GP C ht{g) (resp. hm(( 7 ),hc(g), mdeg((/)) is the head 

term (resp. the head monomial, the head coefficient, the multidegree) of a polynomial g 
so that hm(g) = hc( 5 ) • ht(g) and ht(( 7 ) = hold and ht(GP) = {ht( 5)|5 G GP} 

where R is K,K[t] or K{t). 

Definition 2 (Comprehensive Grobner system (CCS)). Let fix a term order. Let F be 
a subset of (iL[t])[x], Ai ,... ,Ai strata in and GPi ,..., GPi subsets of (iL[t])[x]. A 
finite set Q = {(Ai,GPi),..., (A;,GP/)} of pairs is called a comprehensive Grobner 
system (CGS) on Ai U ■ • • U A; for P if (Ja{GPi) is a Grobner basis of the ideal {(Ja{F)) 
in K[x\ and (ht (era (GP^))) = (ht(GPi)) for each i = 1,... ,l and a G A^. Each (A^, GPi) 
is called a segment of Q. We simply say C/ is a comprehensive Grobner system for F if 
AiU---UA/ = AT™. 


After obtaining a GGS of F w.r.t a total degree term order, as each segment of the GGS 
has the property (ht(cra(GPi))) = (ht(GPi)), the dimension of (GPi) is easily decided 
in Ar[x]. Since an algorithm for computing a CGS terminates, the following algorithm 
clearly terminates. 


Algorithm 1. (ZeroDimension) 

Specification: ZeroDimension(P) 

Testing dimensions of a parametric ideal (P) on P™. 

Input: P: a set of parametric polynomials in (P[t])[x] 

Output: {S, V): S = {(Ai, GPi)..., (A^, GP^)} is a GGS on Ai U • • -UAfe for P s.t. for all 
a £ Ai, (cra(P)) is zero-dimensional in K[x], for each i = I ,.., k.V = {(Afc+i, GPfc+i)..., 
(A;, GP;)} is a GGS on A^+i U ■ • • U A; for P such that for all b G Aj, (a^iF)) is not 
zero-dimensional in K[x], for each j = k + 1,... ,1. = Ai U • • ■ U A*, U A^+i U • • • U A;. 

BEGIN 

5^0; 0;G-(— compute a GGS on P™ of P w.r.t. a total degree term order 

while G 7 ^ 0 do 

select (A, GP) from C; G £- G\{(A, GP)}; d ^compute the dimension of (GP) in K[x] 
if d = 0 then 5 ^ 5 U {(A, GP)} else V <—V U {{A,GP)} end-if 
end-while 
return(5, V) 

END 


In our implementation, we adopt Nabeshima’s algorithm (jNabeshimal . 12012!) for com¬ 
puting comprehensive Grobner systems, because the algorithm is much more useful than 
others for computing dimensions of parametric ideals. 


Definition 3. Using the same notation as in the above algorithm, let (5, V) be an output 
of ZeroDimension(P). Then, for each i = 1... ,k, Ai is called a safe stratum, and for 
each } = fc + 1 ..., Z, Aj is called a danger stratum. 

Example 4. Let / = x\+tx\x\-\-X 2 be a polynomial with a parameter t in (C[t])[xi,X 2 ]. 
A CGS of P = ■§^} w.r.t. the total degree reverse lexicographic term order s.t. 

xi ^ X 2 , is {(V(t),{x?,x|}), (V(t-2), {x\x 2 +xl,x\+xix\}), (V(t + 2), {xfx 2 - 

Xix|}), (C \ V(t(t^ — 4)), {tx\x2 + 2^2, 2xf + tXiX2, — 4)xiX2, — 4 )^ 2 })}. 
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If the parameter t belongs to V(t) or C\V(t(t^ —4)), then (|j, is zero-dimensional. 
If the parameter t belongs to V(t —2) or V(f —2), then §^) is one-dimensional. There¬ 
fore, S = {(V(t), {xf, 3^3 }), (C \ ~4)), {axfx 2 + 2x2, 2x\ + tXiX 2 , — ‘i)xiX 2 , — 

4 ) 3 ^ 2 })} s-iid ^ = {('^(^ ~ 2), {x\x 2 +X 2 , xf + xix\}), {Y{t + 2), {x\x 2 — X 2 ,xf — a:ia;|})}. 
That is, V(t), C\ V(t(t^ — 4)) are safe strata, and Y{t — 2), V(t -I- 2) are danger strata. 

Let (5, X>) denote an output of ZeroDimension(F) where S = {(Ai, GPi)..., 
{Ak,GPk)} and V = {(Afc+i, GPfc-i-i) ■ • ■) (A/, GP/)} (notation is from the algorithm 
ZeroDimension). Since for all a G Ai U • ■ • U A^, {(7a{F)) is zero-dimensional in K[x], 
{aa{F)) is also zero-dimensional in .A[[x]]. However, in general, for all beAk+iU---UAi, 
it is NOT possible for us to say that (i7j(F)) is not zero-dimensional in Ar[[a;]]. For some 
h € Afc+i U • • • U A;, {<Ji{F)) may be zero-dimensional in .A[[a;]]. 

After decomposing the parameter space AT™ into safe strata and danger strata by the 
algorithm ZeroDimension, we compute bases of the vector space Hp on safe strata 
and danger strata, separately. Actually, this decomposition lets us construct an efficient 
algorithm for computing the bases. (See section 3.3). 

As the set F has parameters, the structure of the vector spaces F[f depends on the 
values of parameters t. Here, we introduce a definition of parametric local cohomology 
system oi FIp. 

Definition 5. Using the same notation as in the above, let Aj,Bj strata in A'™ and Si 
a subset of (A'[t]Ai)[^] where I < i < I and I < j < fc. Set S = {(Ai, S'!),..., (A;, S';)} 
and V = {Bi,... ,Bfc}. Then, a pair (<S,T>) is called a parametric local cohomology 
system of on Ai U • • -UAi UBi U • • -UBfc, if for alH G {1,..., 1} and a G Ai, CTo(S'i) is a 
basis of the vector space H„j.(^p), and for all j G {I,..., fc} and h G B^, {c G X\at,{fi){c) = 
■■■ = <^b{fp){c) = 0} ^ {O} where := {V' € k[^] \ (Ta(/i) *-ip = cra(/ 2 ) * V' 

= --- = cralfp)*ip = 0}. 

After here, we represent “a parametric local cohomology system oi FIp on AT™” as 
simply “iJ f” which is the abbreviation. Similarly, we call “a parametric local cohomology 
system oi FIp on a stratum A” “ bases of (the vector space) Flp on A”. 

As this section 3 presents thirteen algorithms for computing bases of the vector space 
Flp, Fig. 1 illustrates the relations of the all algorithms. The main algorithm is ALCo- 
homolog. 

First, we introduce in section 3.2 an algorithm for computing bases of the vector space 
Hp on safe strata. Second, we describe in section 3.3 an algorithm for computing bases 
of the vector space Hp on danger strata. 

3.2. Computation of algebraic local cohomology with parameters on safe strata 

Here, we present an algorithm for computing bases of algebraic local cohomology 
classes Hp, on safe strata. This section consists of three parts. In section 3.2.1, an algo¬ 
rithm for computing monomial elements of bases of iLF is introduced. In section 3.2.2 and 
3.2.3, an algorithm for treating with elements, which form linear combination (^ 
of bases of Hp, is given. 
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Fig. 1. relations of all algorithms 


3.2.1. Monomial elements 

Here, we give an algorithm for computing monomial elements of bases oi Hp. Before 
describing the algorithm, we define some notation. 

Notation 6. Let GP be a set of polynomials in and g € GP. 

(1) The set of monomials of g is denoted by Mono( 5 ), i.e., Mono(g) := {a\x^\g = 

SagN" where a\ G K[t] and a\ ^ 0}. Moreover, the set of monomials of the 
set GP is denoted by Mono(G'P), i.e., Mono(G'P) := Mono( 5 ). 

(2) For all % G {1,... ,n}, a map CV is defined as Changing Variables Xi into ^i. The 
inverse map CV~^ is defined as changing variables into Xi. That is, for any 
g G (iL[t])[a:], CV{g) is in (iF[t])[^]. The set CV{GP) is also defined as CV{GP) = 
{CV{g)\gGGP}. 

For instance, for '^x\x 2 + ^Xi,‘2>x\+Ax 2 G K[xi,X 2 ], then CV(|a;fa ;2 + 5a:i) = |^i ^2 + 
5^1 and CV({|a;fa ;2 + hxi,‘ix\ + 4x2}) = {§^ 1^2 + + 4 ^ 2 } in where 

variables ^i,'f 2 are corresponding to variables xi,X 2 . 

Proposition 7. Let (<S, V) be an output of ZeroDimension(F) and (A, GP) G S. 
Assume that H is a CGS of the monomial ideal (CV(Mono(GP))) in (Ar(t))[^] on A, 
and (A^, G') G B. Then, a monic monomial tp = ^“’^■^2 ^ ''' Cn" which does not belong to 
(ht(G')), has the property fi^ip = 0 for each 1 < z < p. Namely, all terms which do not 
belong to (ht(G')), are members of bases of Hp on A^ 

Proof. Let ip = ''' Cn" be a monomial s.t. ip ^ {ht{G')). By Definition 2, for 

all a G A', (ht(G')) = {ht((Ta(G'))) = (CV(Mono( cra(GP)))). As (ht(G')) is a zero¬ 
dimensional ideal and ip ^ {ht(G')), for all G CV( Mono(C)), there always 

exists j G {1, 2,..., n} such that Xj > aj where (oi, 02 , • ■ •,««), (Ai, A 2 ,..., A„) G N". 
Therefore, by the multiplication, fi * ip = 0. □ 









































































This proposition gives rise to the following algorithm to compute monomial elements 
of bases oi Hp on A. Since the termination of Nabeshima’s algorithm (INabeshimal . l2012r) 
is guaranteed, the following algorithm terminates. 


Algorithm 2. (MonoSafe) 

Specification: MonoSafe(A, GP) 

Computing monomial elements of bases of on a safe stratum A. 

Input: (A, GP): a segment of a CCS of F such that for all a € A, {aa{F)) is zero¬ 
dimensional in K\x\. (This is from ZeroDimension(P).) 

Output: Ad : a finite set of triples (A^, M, G) such that the set M includes all monomial 
elements of bases oi Hp on A', and the elements of M do not belong to (G). 

BEGIN 

Ad 0 ; B ^ compute a CGS of CV(Mono(GP)) on A 

while P 7 ^ 0 do 

select (A', G') from B- B ^ B\{{M , G')}; G ^ ht(G') 

M •<— compute monomial elements which do not belong to (G) in K[^ (=i=l) 

Ad ^AdU{(A',M,G)} 

end-while 
return Ad 
END 


Let us remark that as (GP) is a zero-dimensional ideal on A, the set M consists of 
finitely many monomial elements. Note that monomial elements, on danger strata, will 
be considered in section 3 . 3 . 

We illustrate the algorithm MonoSafe with the following example. 

Example 8. Let / = xf+tx\x2 + X2 be a polynomial with a parameter t in (C[t])[a;i,a:2]. 

^ Then, F satisfies generically {a € X\^{a) = ^{a) = 0 } = {O} 

where A is a neighborhood of the origin O of C^. From Example 4 , (V(t), {xi,X2}) and 
(C \ V(t(t^ — 4 )), {txlx2 + 2 ^ 2 , 2 xi+txiX 2 j {t^ —A)xiX 2, {t^ —4)3:2}) can be inputs of the 
algorithm MonoSafe. 

(1) Take (V(t), {xf, xf}) as an input of the algorithm MonoSafe. Then, a CGS of 

CV({x?,xi}) on V(t), is {(V(t),U?,e|})}. Set G, = ht({e?,e|}) = Then, 

all elements of Mi = { 1 ,Ci,6,Ci,66,^2.C16,6C2:do not belong to (Gi). 
See • in Fig. 2. Mi can be a subset of bases of Hp on V(t). 

(2) Take (C\V(t(t^— 4 )),P) where GP = {tx\x2 + 2 x 2 , 2 x\+txix\, — A)xiX2, — 

4)x2}. As Mono(GP) = {txiX2, 2 x 1 , 2 xi, teix^, (t^ — 4)xiX2,(t^ —4)x2}, a GGS 
of (CV(Mono(GP))) on C\ V(t(6 - 4 )) is {C\V(t(t 2 - 4)), {el^^66,^66^6"})}• 
Set G2 = ht({6,t66,^66,C2}) = {6,Ci6, 66,61 and compute monomial 
elements which do not belong to (G2). Then, we obtain M2 = { 1 , 6 , 6 , 6 , 66 , 

which can be a subset of bases of F[p on C\V(<(t^ — 4 )). See • in Fig. 3 . 
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3.2.2. Head terms of linear combination elements and the main algorithm 

Here, we illustrate an algorithm for computing bases of Hp- Before describing the 
algorithm, first we treat with elements, which form linear combination (^ cxf^), of bases 
oiHp. Especially, we discuss how to decide head terms of the linear combination elements 
(^ c\f^). Second, an algorithm for computing bases of Hp on safe strata, is given. Note 
that an algorithm for deciding lower terms, will be described in section 3.2.3. 

Let us recall the following lemma which follows from the fact that if '0 G Hp, so is 
Xi * tp € Hp for each i = 1,..., n. This lemma informs us candidates of head terms in 
Hp. 


Lemma 9 ( Taiima and Nakamura ( 2009ll l. Let Ap denote the set of exponents of head 
terms in Hp and A = (Ai,..., A„) G N". Let denote a subset of Ap : Ap = {A G 
N" Idi/) G Hp s.t.ht('0) = and = {A' G Ap\X' -< A}. If A G Ap, then, for each 
j = 1,2,... ,n, {Xi, X 2 ,..., Xj-i, Xj - l,Xj+i,...,Xn) is in A^^^ provided Xj > 1. 


Let be a term where A = (Ai,..., A^) G N”. We call ■ fi a neighbor of for 
each i = l,...,n. Then, • fi\ = + 1- 


Notation 10. Let T be a set of terms in K\^. Then, we define the neighbor of T as 
Neighbor(r), i.e., Neighbor(r) := {r • ^i|r G T, i = 1,..., n}. 


The following corollary is a direct consequence of Lemma 9. 

Corollary 11. Let TList*-'^^ = {^''|A G Af,!^'''] = d}. If for all i G {!,..., n}, t = 
■■■'?«" G TList^'^'*'^^ satishes ^i|r, then, {r/^ijAi 0,i G {l,...,n}} C TList^'^^ 
where (Ai,..., A„) G N". 

Let T be a subset of TList^'^^ Then, by Corollary 11, there is a possibility that an 
element of Neighbor(T) belongs to TList^‘^~'’^^ This fact makes up the following algo¬ 
rithm which outputs new candidates for head terms. 


Algorithm 3. (HLem) 

Specification: HLem(T, TList*-'^^) 

Making new candidates for head terms from T. 

Input: T: a set of terms whose total degree are d, and T C TList^'^^ 
Output: S: a set of new candidates whose total degree are d -I- 1. 
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BEGIN 

S-5— 0; B<— Neighbor(T) 
while B 7 ^ 0 do 

select r from B; B-(— B\{t} 
for i from 1 to n do Flag-S— 1 
if ^i|r then 

if r/^i ^ TList^*^^ then Flag'(— 0 ; break end-if 
end-if 
end-for 

if Flag= 1 then S'(— SU{r} end-if 
end-while 
return S 

END _ 

If r is not in the set of head terms oi Hp (written ht(iJF)), then neighbors of r are 
not in ]it{Hp). This means that if r is not in htlHp), then it is unnecessary to compute 
elements which are divided by r, as candidates for head terms. This fact makes up the 
following notation NonMember. We also give the notation Car and Cdr which are 
exploited in the some algorithms. 

Notation 12. Let T be a set of terms in and r be the smallest element in T w.r.t. 
the term order (of Definition 1). 

(1) Let FL be a set of terms in such that for all € FL, A is not in where 
Af is the set of exponents of head terms in Then, the notation NonMember 
of (T, FL) is defined by NonMember(r, FL) := {tp € T\ip \ tp for all (p € FL}. 

(2) The notation Car and Cdr for T, are defined as follows 

Car(T) := € T I IV'I = |t|} , Cdr(r) := T\Car(T). 

(3) Suppose that e T\\^^\ = d € N} and TT = ..., 

where € N and ^ 0 for each i = 1,... ,u. Let dj be the minimal number in 
{di,..., d„}. Then, the same notation Car and Cdr are defined, for a set of sets 

TT, as follows ,. , 

Car(TT) := Cdr(TT) := TT \Car(TT). 

Let A4 be an output of MonoSafe(A, GF) where (A, GF) is a segment of a CCS of F. 
Suppose that (A', M', G') is an element of M.. Remember that all elements of M' do not 
belong to {G'). Since clearly G' C Neighbor(M'), elements of G' become candidates 
of head terms in Hp. The use of this property makes candidates of the head terms, 
efficiently. 

Corollary 13. Using the same notation as in the above discussion. Notation 12 and 
Lemma 9, let € M', = dj and = TList*-‘^^\M*^‘^). Then, elements 

of Neighbor(T*^‘^)) and G' can be candidates of head terms in Hp. 

Suppose that GList = ..., and FL is a set of terms in such that 

for all G FL, A is not in Af where G G'||^'>'| = d^j and j G {!,.■. ,u}. 

Now, we introduce how to obtain a set of candidates of head terms in Hp from and 
GList. In order to make the set of the candidates, the following four cases are considered. 


11 




Case (i) = 0 A GList = 0. Case (ii) = 0 A GList + 0. 

Case (iii) ^ 0 A GList = 0. Case (iv) 0 A GList ^ 0. 

In case (i), our main algorithm terminates, because any candidates of the head terms 
can not be made by the sets. In case (ii), Car (GList) has to be considered as a set of the 
next candidates w.r.t. the term order. In case (iii), NonMember(HLem(r(‘^\ TList^'^^), 
FL) has to be considered as a set of the next candidates whose total degree is fi+1. In case 
(iv), for any r GCar(GList), if |t| — d = I, then NonMember(HLem(T*^‘^\ TList*-"^^), 
FL)UCar(GList) has to be considered as a set of the next candidates, otherwise the next 
candidates is NonMember(HLem(r(‘^\ TList*-"^^), FL). 

Let us remark that the algorithm BodySafe decides head terms of bases of Up, from 
bottom to up with respect to the term order (total degree lexicographic term order). 
Therefore, the sets and TList*-*^^ are already obtained when the following algorithm 
makes the set of the candidates whose total degree are d + I. 

Algorithm 4. (HeadCand) 

Specification: HeadCand(T(‘^\ GList, TList*^‘^\ FL) 

Making new candidates for head terms. 

Input: GList, TList^-'^^ FL: described above. 

Output: CT: a set of new candidates for head terms (or Car(GList)). GList: renewed 
GList; r('*):renewed 

BEGIN 

if 2i(d) _ 0 A GList = 0 then /* case (i) */ 

GT ^ 0 ; return(CT, GList, T^'^)) 
else if = 0 A GList 0 then /* case (ii) */ 

GT ^ Car(GList); GList Cdr(GList) ; return(GT, GList, 
else if 7 ^ 0 A GList = 0 then /* case (iii) */ 

GT ^ NonMember(HLem(r(^), TList^*^^), FL); d ^ d + 1; ^ 0 

return(CT, GList, r(^)) 

else if 7 ^ 0 A GList 7 ^ 0 then /* case (iv) */ 

G <—Car(GList) ; GL'(— Cdr(GList); select from G 

if — d > I then 

GT ^ NonMember(HLem(r(^), TList^'^^), FL); d ^ d + 1 ; ^ 0 

return(CT, GList, r(^)) 

end-if 

if — d = I then 

GT ^NonMember(HLem(T(‘^),TList(‘')), FL)UG ; GList ^GL 
return(CT, GList, r(^)) 

end-if 

end-if 

END _ 

The algorithm BodySafe consists of mainly two parts, computing candidates for head 
terms and lower terms. For each part, the algorithm makes use of sets as intermediate 
data. As this is a dynamic algorithm, each intermediate data is often renewed in the algo¬ 
rithm. As sets SList, MList, LList, GList, GT, GL are frequently used in algorithms 
on a stratum, we fix the meaning of the sets as follows. 
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Notation 14. SList := {'ijj € is a linear combination element of a basis}. 

MList := {tp G K[^]\ ip is a monic monomial element of a basis}. 

LList := G is a lower term of ip where ip G SList}. 

CT := {t G T is a candidate for head terms of a basis}. 

FL := {t G T is a failed candidate for head terms}. 

GList := lJj{{C^ G G\ = d^}} described in the algorithm HeadCand. 

:= {r G r is a head term whose total degree is d}. 

CL := G is a candidate for lower terms for some t G CT}. 

As F has parameters, when we compute bases of by the main algorithm ALCoho- 
mology, the parameter space A™ is decomposed to suitable strata for the bases. Hence, 
on each stratum, the sets above are decided. Note that when the algorithm terminates, 
then a set SList U MList becomes a basis oi Hp on each stratum. 

In the following two algorithms, sets EL,LL,UU, RR are used for algorithmic consis¬ 
tency, to decide lower terms. The sets will be explained in section 3.2.3. 

The main algorithm ALCohomology consists of two parts for safe strata and danger 
strata. The first part an algorithm BodySafe for safe strata, is given in this section. The 
second part an algorithm BodyDanger for danger strata will be discussed in section 3.3. 

Suppose that Q is a list. Then, Q[i] means the ith element of the list Q. For example, 
let Q = [A, CT, GList], then Q[l] = A, Q[2] = CT and Q[3] = GList. In the following 
algorithms, lists Q, £ and MList^'^^ = {r G MList ||t| = d} play actively. 


Algorithm 5. (ALCohomology) 

Specification: ALCohomology(F, k) 

Computing bases of a vector space Hp with parameters. 

Input: F = {/i,..., /p}: F C {K[t])[x\ satisfying generically {a G Aj /i(a) = • • • = 
fp{a) = 0} = {0} where A is a neighborhood of the origin O of A", G N: an estimated 
bound of dimensions of the vector space Hp or a sufficient big number (see section 3.3). 
Output: (5, X>): 5 is a set of lists [A, SList, MList, LList, FL] where SList U MList is a 
basis of Hp on A, LList is a set of lower terms of SList and FL is a set of failed candidates 
for head terms on A. 

BEGIN 

CT ^ 0; SList ^ 0; LList ^ 0; FL ^ 0; LL ^0; RR ^ 0; EL ^ 0 
UU ■<— 0; AC G- 0; VC g- 0; {Z,M) G- ZeroDimension(F) 

/*on safe strata */ 
whileZ ^ 0 do 

select Zi from Z G- Z\{Zi}; A4 G- MonoSafe(Zi) 
while AI 0 do 

select (A, MList, G) from M.; M g~ AI\{(A, MList, G)} 

GList ■(— S = di}}; t g- the smallest element in G w.r.t. -<; d G- |t| 

TA) ^ 0; g ^ [A, CT, GList, SList, MList, LList, FL, LL, EL, RR, UU] 
AC^ACU{Q} 

end-while 

end-while 

Coho G~ BodySafe(AC, F) 
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/*on danger strata */ 
whileAA 7 ^ 0 do 

select Ni from A/”; J\f ^ J\f\{Ni}] A4 -s— MonoDanger(A^i) 
while Ad 0 do 

select (A, MList, G) from Ad; Ad ^ Ad\{(A, MList, G)} 
if MList 7 ^ 0 then 

GList ■(— ^ I = di}}; t ^ the smallest element in G w.r.t. -<; d ■(— |t| 

TG) ^ 0. g ^ [A, CT, GList, T^, SList, MList, LList, FL, LL, EL, RR, UU] 

VC ^ VCiJ{Q] 
else V i— VU {A} 

end-if 

end-while 

end-while 

(Co,X>i) <— BodyDangerl:^,X>£,F) 

5 ^ CohoU Go; V^VUVi 
return (S, V) 

END 


Algorithm 6. (BodySafe) 

Specification: BodySafe(AC,F) 

Computing bases of algebraic local cohomology Hp for AC. 

Input: AC: a set of lists ([A, CT, GList, SList, MList, LList, FL, LL, EL, RR, UU]). 
Output: S: a set of lists [A, SList, MList, LList, EL] where SList U MList is a basis of Hp 
on A , LList is a set of lower terms of SList, and FL is a set of failed candidates for head 
terms on A. 

BEGIN 

5e-0 

while AC 7 ^ 0 do 

select £ = [A, CT, GList, SList, MList, LList, FL, LL, EL, RR, UU] from AC 
AC ^ AC\{£:} 

if CT ^ 0 then e- Car(CT); CT ^ Cdr(CT) 
else 

(CT, GList, c-HeadCand(TW, GList, MList^'^^UT^''), EL) (^1) 
if CT 7 ^ 0 then 

^ Car(CT); CT ^ Cdr(CT) 
else 

S c- 5 U {[A, SList, MList, LList, EL]} 

end-if 

end-if 

(CL, UU, EL) 4- LowCand(^T', SList, MList, LList, LL, UU, RR, EL) (^2) 

Q ^ [CT, GList, MList, UU] 

V ^ OneElement(gT', CL, A, EL, EL, SList, LList, Q, F) _ (03) 

iS <S U BodySafe(F) 

end-while 

return S 

END 
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The algorithm BodySafe consists of three parts (<C>1), (<C>2) and (<C>3). In (<C>1), new 
candidates for head terms are computed. The part (<C>1) was already described in the 
beginning of this section. In (<C>2), candidates (CL) of ^’’'’s lower terms are computed. 
The part (<C>2) will be described in section 3.3. Here, we do not explain the part (<C>2), but 
by seeing the operation of CL, one can understand the flow of the algorithm BodySafe. 
In (03), an element + X^AeCL^ ^(^)) is tested whether it can be in Hp 
or not. That is, linear combination elements are decided in the part (03). Note that in 
(03), a list Q is not essentially used by the algorithm OneElement. The list Q is just 
used in order to shorten the algorithm. The part (03) is given as follows. 


Algorithm 7. OneElement 

Specification: OneElement(^^, CL, A, T(^), EL, FL, SList, LList, Q, {/i,..., fp}) 
Testing whether + X^AeCL is in Hp or not. 

Input: CL, A, EL, FL, SList, LList, Q: described in the algorithm BodySafe. 

Output: C: a set of lists [A, CT, GList, T(^), SList, MList, LList, FL, LL, EL, RR, UU]. 
BEGIN 


0 ^ set r +EeGCL caC^ where ca’s are indeterminates 

for i from 1 to p do 

i’ /*check /i * 0 = 0. /i * 0 G {K[t, ca])[C]*/ 

while 0 7 ^ 0 do 

E E Li {hc(0) = 0}; — hm(0) 

end-while 

end-for 


(Ai, A 2 ) solve the system E of parametric linear equations on A. 


(* 1 ) 


while Ai 7 ^: 0 do 

select an element (A', [ca’s solutions]) from Ai; Ai -s— Ai\{(A', [ca’s solutions])} 
f/'' -5— substitute ca’s solutions into - ip ] SList -S— SList U{'i/''}; TA) tA) u 
(EL, LL, RR, LList) ^ renew_low(l, EL, ip' — , LList) 

£ 7- £ U {[A', Q[l], Q[2], TA), SList, Q[3], LList, FL, LL, EL, RR, Q[4]]} 


while-end 
while A 2 7 ^ 0 do 

select an element A'from A 2 ; A 2 A 2 \{A'}; FLFL U{^'''} 

(EL, LL, RR, LList) renew_low(0, EL, LList) 

£ <- £ U {[A', Q[l], Q[2],tA), SList, Q[3], LList, FL, LL, EL, RR, Q[4]]} 


while-end 
return £ 


END _ 

If V' = ?^+Ee6CL caC^ is in E[p, then ip satisfies conditions fi * ip = 0 for each 
i = I,.. . ,p. These conditions give us a set E of ca’s linear equations. Thus, by solving 
the system E, we know whether ip is in Elp or not. Namely, if solutions of ca’s exist, then 
Ip G Hp, and if the solutions of ca’s do not exist, then ip ^ Hp. 

Let us remark that as the system of equations E has parameters, the stratum A has 
to be decomposed into suitable strata for the solutions. For instance, let t be a parameter 
and X, y be variables. Consider a system “tx + y = 4, 3x + 2y = —9” of parametric linear 
equations on C\V(t). Then, the system has the following solutions; if the parameter 
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t belongs to C\V(t(3< + 4)(2t — 3)), then x = 2 ^,y = 2 t- 3 ^ ’ 


if the parameter t 


belongs to V(3t + 4), then x = —3, j/ = 0, and if a parameter t belongs to V(2t — 3), 
then E has no so lution. There exist several algori thms for solving a system of parametric 
linear equations (jCao and Choi]l . Il992l: ISitl . Il992l l. In our implementation, we extend the 
Gaussian elimination method to handle parametric cases. 

In the box (^I) of the algorithm OneElement, Ai means a set of pairs (A^, [ca’s 
solutions]) and A 2 means a set of strata such that for any stratum of A 2 , the system 
has no solution. The algorithm OneElement has a subalgorithm renewJow which 
computes candidates of lower terms of and is given in section 3.2.3. 


Theorem 15. The first part of the algorithm ALCohomology (i.e., BodySafe) ter¬ 
minates and outputs a set Coho which has a list [A, SList, MList, LList,FL] such that 
SList U MList is a basis oi Hp on A. 


Proof. The algorithms LowCand and renewJow are considered in section 3.2.3. and 
the termination and correctness are discussed in section 3.2.3. In the algorithm Ze- 
roDimension, the parameter space A™ is decomposed to a finite number of strata. 
As we described, in the algorithm OneElement, an al gorithm for solv i ng the sy s tem o f 
parametric equations, outputs a finite number of strata (jCao and ChoJ . ll992l : ISitl . 1199211 . 
Since the algorithm BodySafe works on safe strata, (F) is zero-dimensional on the s trata . 
This means that H p is a finite-dimensional vector space ( Taiima and Nakamural 2009t 
iTaiima et al.l . I 2 OO 9 II . Therefore, the first part of the algorithm ALCohomology (i.e., 
BodySafe) generates a hnite number of strata. Thus, the algorithm terminates. More¬ 
over, clearly all elements of SList U MList are linearly independent on A, SList U MList 
is a basis of Hp on A. □ 


Example 16. Let f = xf + tx^x"^ + X 2 G (C[t])[a::i, 3 : 2 ]. Set F = {■§^, Then, F 
satisfies generically {a € A|^-(a) = ■§^{a) = 0} = {O} where A is a neighborhood of 
the origin O of . The term order is the total degree lexicographic term order such that 
6 -< 6 - 

(0) CT = 0, SList = 0, FL = 0. 

Case (1): From Example 8 , the output of MonoSafe(V(t), {ccf, a^f}) is (V(t), MList, Gi) 
where MList = {1,Ci, 6 ,Ci, 66 ,^ 2 .C 16 , 6 C 2 :C 1 C 2 } and Gi = { 6 , ^f}. Then, GList = 
{{Cl, 6 }} and the smallest element in Gi is 6 - Since | 6 l = 3, set = 0. 

(1-1) As = 0, GList = {{Cl,Cl}} and CT = 0, CT and GList are renewed as 
CT = Car(GList) = { 6 , 6 } and GList = 0. So, Car(CT) = 6 and CT = 
Cdr(CT) = { 6 }- Then, the algorithm LowCand outputs the empty set as the set 
of candidates of 6 ’s lower terms. This means that there is no candidate for lower 
terms of 6 - Since *6=4, ■§^ *6 = 0, we obtain equations 4 = 0 and 0 = 0. 
Clearly, 4 = 0 is false. Hence, 6 can not be a member of Hp. Renew FL as { 6 }- 
(1-2) As CT = { 6 }, the new candidate for head terms is 6 - Renew CT as Cdr(CT) = 0. 
Then, the algorithm LowCand outputs { 6 } as the set of candidates for lower 
terms. Set 6 +C( 3 ,o)Ci where C( 3 ^o) is an indeterminate. Since * (6 + C( 3 ,o)Ci) = 
4 c( 3 ,o), * (Cl + C( 3 ,o)Ci) = 4, we obtain equations 4 c( 3 _o) = 0 and 4 = 0. Clearly, 

4 = 0 is false. Hence, 6 can not be a head term in Hp. Renew FL as { 6 , C|}- This 
process terminates, because CT = GList = 0. 
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Case (2): From Example 8, the output of MonoSafe(C\V(t(t^— 4 )), {tx\x2 +2x2, 2 xf + 

txixi, 4 )xiX2, 4)x 2}), is (C\ 4 )), MList, G2) where MList = { 1 , Ci, $2, 

C16, ^ 1 } and G2 = Ui,??6,?iC|,? 2}- Then, GList = and the 

smallest element in G2 is Since = 3 , set = 0 . 

( 2 - 1 ) As = 0 , GList ^ 0 and GT = 0 , GT and GList are renewed as CT = 
Car(GList) = {C?, 6 ^ 1 , ^ 1 } and GList = 0 . So, Car(GT) = and CT = 

Cdr(CT) = {^1^2,CiCl)Cl}- By the same reasoning as in ( 1 - 1 ), Ci can not be a 
member oi Hp. Renew FL as {Cil- 

( 2 - 2 ) As CT = {CiCli C1C2, ClI) the new candidate of a head term in Hp, is CiC2- Renew 
CT as Cdr(CT) = {CliCiCl}- Then, the algorithm LowCand outputs {Ci} as 
the set of candidates of CiC2’s lower terms. Set C1C2 + C(3,o)Ci where C(3^o) is an 
indeterminate. Since * (C?C2 + C(3,o)Ci) = 4 c( 3 ,o), ^ * (C?C2 + C(3,oCi) = 2 t, a 
system of equations is “ 4 c( 3 ^o) = 0 , 2 t = 0 ”. As we work on the stratum C \ — 

4 )), 2 t = 0 is false. Hence, C1C2 can not be a head term in Hp. FL = {Ci}U{CiC2} = 

{Cl,CfC2}. 

( 2 - 3 ) As CT = {CiCIjCI}) the next candidate is CiCl- Renew CT as Cdr(CT) = {ClI- 
Then, the algorithm LowCand outputs {Ci, C1C2} as the set of candidates of CiCl’s 
lower terms. Set ^ = CiC2+C(2,i)CiC2+C(3,o)Ci where C(2p), C(3^o) are indeterminates. 
Since = 4c(3_o) + 2 t, ■§^*’ 4 ’ = 2tc(2,i), a system of equations is “4c(3_o) + 2 t = 

0 , 2te(2p) = 0 ”. Solve the linear equations, then C(3^o) = and C(2,i) = 0 . Hence, 
C1C2 “ jtCi is a member of a basis of Hp on C\V(t(t^ — 4 )). Thus, SList = 
{CiCl - = {C1C2}. 

( 2 - 4 ) As CT = {ClI, the next candidate is Cl- Renew CT as Cdr(CT) = 0 . Then, 
the algorithm LowCand outputs {C17C1C2} as the set of candidates of Cl’s lower 
terms. Set V' = Cl + C(2,i)CiC2 +C(3,o)Ci where C(2,i), C(3_o) are indeterminates. Since 
^ = 4 c( 3 ,o), ^ * V' = 2 fc( 2 ,i) + 4 , a system of equations is “ 4 c( 3 ,o) = 0 , 

2 tc( 2 ,i) +4 = 0 ”. Solve the linear equations, then C(3_o) = 0 and C(^2,i) = ~^- 
Hence, C| — ^CiC2 is a member of a basis of Hp on C \ V(t(<^ — 4 )). Thus, SList = 
{CiCl - ltd, cl - ^dC2}, T( 3 ) = {C1C2", c|}- Now, as CT = 0 , the set CT should be 
renewed by the algorithm HeadCand. Since Neighbor(T^^)) = {C1C2,CiCl,ClI 
and C1C2IC1C2 where C1C2 is in FL, C1C2 can not be in CT. Therefore, the renewed 
CT is {CiCl.ClI- 

(2-5) As CT = {CiCliCll, the next candidate is CiCl- Renew CT as Cdr(CT) = {Cll- 
Then, the algorithm LowCand outputs {Ci, C 1 C 2 , C 1 C 2 , CiCliCil as the set of can¬ 
didates of CiCl’s lower terms. Set V’ = C 1 C 2 + C( 2 , 2 )CiC 2 + C( 3 .i)CiC 2 + C( 4 ,o)Ci + 
C(2 ,i)CiC2 +C(3,0 )Ci wherec(2,2),C(3,i),C(4,0),C(2,i),C(3,0) are indeterminates. = 

(4C(4_o)+ 2te(2,2))Cl + (4C(3 p)+ 2t)C2+4C(3 0 ), = (2tC(3 i)+4)Cl+2tC(2_2)C2 + 

2tc(2,i). If V’ is in Hp, then ip satisfies the conditions *p = 0 and *p = 0. 
Hence, we have to check the five equations 4c(4 g) + 2te(2,2) = 0, 4 c(3 1 ) + 2t = 
0, 2fc(3^i) + 4 = 0, 2tc^2,2) = 0, 2te(2p) = 0 on c'\V(t(f^ ~4)). The two equations 
4 c(3 1 ) + 2t = 0 and 2fc(3 1 ) + 4 = 0 hold only if t = ±2. Therefore, p) can not be in 

Hp'. FL = {ci,dC2} u {dcii = {c?,dC2,CiC|}- 

( 2 - 6 ) The next candidate is Cl and CT = 0 . The algorithm LowCand outputs {Ci,CiC2, 
Cl,C1C2,C1C2,CiCll as the set of candidates of Cl’s lower terms. Set = Cl + 
C(l,3)ClCl + C(2,2)ClC2 + C(3,1 )CiC 2 + C(4,o)Cl + C(2,l)ClC2 + C(3,o)Cl where C(1_3), C(3_l), 
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C( 4 ,o).C( 2 ,i):C( 3 ,o) are indeterminates. * ip = ( 2 fc(i_ 3 ) + 2 te( 2 , 2 ) + 4 c( 4 _o))^i + 
4 c( 3 _i )^2 + 4 c( 3 _o)) = ( 4 c(i_ 3 ) + 2 te (3 1 ))^ + ( 2 te( 2 _ 2 ) + 4)^2 + 2 ic( 2 ,i)- Hence, 

we have to check the system of equations: 2 tc(i_ 3 ) + 2 te( 2 , 2 ) + 4 c( 4 _o) = 0 ; 4 c( 3 _i) = 
Oj 4 c( 3 _o) = 0, 4 c(i_ 3 ) + 2 te( 3 ^i) = 0, 2tc(2,2) +4 = 0, 2 te( 2 q) = 0. Then, the 
solution is : C(i_ 3 ) = 0 ,C( 2 , 2 ) = - 7 ,C( 3 .i) = 0 ,C( 4 _o) = l,C( 2 .i) = 0 ,C( 3 _o) = 0 . 
Hence, ~ + Ci is a member of a basis of Hp on C\V(t(t^ — 4)). Thus, 

SList = {^ 1^1 - it??,el - - i???i + ??}, rW = {^ 4 |_ As CT = 0, CT 

can be renewed as {ef}- 

(2-7) The next candidate is ?! and CT = 0. The algorithm LowCand outputs {?i,?i? 2 , 
?i j ?i? 2 , ?i? 2 ) ?i? 2 ) ?i} as the set of candidates of ef’s lower terms. Set V' = ?! + 
C( 5 ,o)?i + C(i, 3 )?i?| + C( 2 , 2 )?i ?2 + C( 3 .i)?i ?2 + C( 4 ,o)?i + C( 2 .i)?i ?2 + C( 3 _o)?i where 
C( 5 ,o)i C(i, 3 )! C( 2 , 2 )i C( 3 ,i)i C( 4 ,o)! C( 2 ^i)C (3 o) are indeterminates. In this case, there is 
no solution that satisfies the conditions if ip = 0 and * ip = 0. FL = 
{?!,?i? 2 ,?i? 2 )el}- t^s CT = GList = 0, this process terminates. 

We summarize the results as follows: 

if a parameter t belongs to V(f) (i.e., t = 0), {1,?!,e2,?!,eie2,?!)eie2,ei?!:eieil is 
a basis ot Hp (algebraic local cohomology classes), 
if a parameter t belongs to C\ V(t(t^—4)) (i.e., t ^ 0, t 7 ^ ±2), then {1, ?!; e 2 j ?!! eie 2 , 

el,?!?! - ^^ei,?! - ^e?e2,el - ie??!+e?} is a basis of Hp. 

In Fig. 4 and 5, we represent an element of MList as • and an element of ht(SList) as 
A. Note that on each stratum, a basis of Hp is MList U SList. As the set FL plays a key 
role to construct standard bases (see section 4), we specially give the elements of FL in 
the Figures. 




3.2.3. Lower terms of linear combination elements 

The aim of this section is to construct subalgorithms “LowCand” and “renewjow” 
which are in the algorithms “BodySafe” and “OneElement” . Here, we discuss how to 
compute candidates of lower terms. The ideal for computing the candidates efficiently, is 
to use the information of the intermediate data SList, MList, LList, FL. Before describing 
the algorithms, we introduce the following useful lemma. 


Lemma 17 ( Taiima and Nakamural ( 2009ll '). Using the same notation as in Lemma 9, 
let Ap denote the set of exponents of lower terms in Uf and denote a subset of 
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Ap : = {A' G Ai?|A' ^ A}. If A G Ap, then, for each j = 1,2,..., n, (Ai, A 2 ,..., 

Aj_i, Aj — 1, Aj+i,..., A„) is in A^^ U provided Xj > 1. 

The algorithm BodySafe computes linear combination elements of a basis of Hp from 
bottom to up with respect to the term order. The next corollary shows a relation between 
the indeterminate data “SList, MList, LList” and new candidates of lower terms. 

Corollary 18. Let A = (Ai,..., A„) G N" and let SList, MList, LList be indeterminate 
data in the algorithm BodySafe. If G LList, then, for each j = l,2,...,n, the 
term is in ht(SList) U MList U LList, provided Xj > 1. Conversely, an element of 

Neighbor(ht(SList) U MList U LList) becomes a candidate for lower terms. 

This corollary leads us to construct the following algorithm which is essentially same 
as the algorithm HLem. 

Algorithm 8. (LLem) 

Specification: LLem(Ne, SList, MList, LList) 

Making candidates for lower terms from Ne. 

Input:Ne: a set of terms. 

Output: S: a set of new candidates for lower terms. 

BEGIN 

S'^0 

while Ne ^ 0 do 

select T from Ne; Ne Ne\{T} 
for i from 1 to n do Flag^ 1 
if ^i|r then 

if r/^i ^ ht(SList) U MList U LList then Flag'(— 0; break 
end-if 
end-if 
end-for 

if Flag= 1 then S' ^ S' U {t} end-if 
end-while 
return S 

END _ 

Let us remark that if a lower term is in ht(SList) U MList, then the lower term can be 
reduced by elements of SList U MList. Namely, LList obtained, becomes always a part of 
candidates for lower terms. Thus, a set of the candidates is 

CL = {proper new candidates of lower terms} U LList. 

As sets EL, LL, RR, UU are frequently used in the algorithms LowCand and re- 
newJow on a stratum, we fix the meaning of the sets as follows. 

Notation 19. EL := G R1[C]| is a new candidate for lower terms. ^ LList}. 
LL := is a proper new lower term which belong to EL}. RR := EL\LL. 

UU := G Neighbor(LL)| ^ for some G CT}. 

Note that a set EL U LList becomes a set of candidates for lower terms. 


19 






As LL is a set of new lower terms, by Corollary 18, elements of Neighbor(LL) become 
candidates of lower terms. Furthermore, elements of ^ G UU} also become 

candidates of lower terms where is a candidate of a head term. 

Algorithm 9. (LowCand) 

Specification: LowCand(^''', SList, MList, LList, LL, UU, RR, EL) 

Making candidates for lower terms of 

Input: SList, MList, LList, LL, UU, RR, EL: described in BodySafe. 

Output: (CL,UU,EL): elements of CL are candidates for ^'’'’s lower terms. UU is a 
renewed set. EL is a renewed set. 

BEGIN 

guu} 

if LL = 0 then 

LU ^ LLem([/, SList, MList, LList); UU ^ UU \U 
EL ^ EL \LU; CL ^ LList U EL 
return(CL, UU, EL) 

else /* make EL, UU from RR, LL .*/ 

UU ^ (UU\f7)\{^T'}; RR^UURR 
B ^ {/3| C'*' ^ /3 G Neighbor(LL))}; UU ^ B U UU 

D ■(— LLem(Neighbor(LL)\i?, SList, MList, LList) 

EL^ (£i\(i:)nRR))URR; CL ^ ELULList 
return(CL, UU, EL) 
end-if 
END 

Since the algorithm OneElement is dynamic, each intermediate data of EL, LL, RR 
and LList, is often renewed in the algorithm. If a system of linear equations has solutions 
(i.e., Z = 1 in renewJow), then the proper new lower terms appear as LL. If a system of 
linear equations does not have any solution (i.e., Z = 0 in renewjow), then the candi¬ 
date of a head term becomes a candidate of lower terms because the candidate is always in 
Neighbor(ht(SList)UMList U LList). This observation makes the algorithm renewJow. 

Algorithm 10. (renewjow) 

Specification: renew Jow(Z, EL, ip, LList) 

Renewing the sets EL, LL and RR. 

Input: Z: 0 or I. ip: a polynomial. 

Output: (EL, LL, RR, LList): renewed sets EL, LL, RR, LList. 

BEGIN 

if Z = 0 then LL <— 0; EL ELU{7/’}; RR 0 
else LL ^ ht(Mono(^/’)) fl EL; LList LList U LL 
if LL 7^0 then RR^EL\LL; EL ^ 0 

end-if 

end-if 

return (EL, LL, RR, LList) 

END 
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Example 20. Let us consider Example 16, again. Here, we show the process for comput¬ 
ing candidates for lower terms according to the algorithms LowCand and renewjow. 
Since a set of candidates for lower terms is CL, we mainly observe the set CL. 

(0) LList = 0, LL = 0, RR = 0, EL = 0, UU = 0. 

Case (1): First, we start to discuss how to compute CL on V(t). 

(1-1) Take as a candidate of a head term. The algorithm LowCand outputs the empty 
set as CL. By the algorithm renewJow, the set EL is renewed as {^i}. 

(1-2) Take a candidate of a head term. According to the algorithm LowCand, 

CL = EL U LList = {^f}. In Example 16, can not be a head term in Hp. Hence, 
EL is renewed as 

Case (2): Second, we discuss how to compute the set CL on \ V(t(t^ — 4)). 

(2-1) Take as a candidate of a head term. The algorithm LowCand outputs the empty 
set as CL. By the algorithm renewJow, the set EL is renewed as {^i}. 

(2-2) Take ^^^2 as a candidate of a head term. According to the algorithm LowCand, 
CL = EL U LList = {■Ci}- In Example 16, can not be a head term in Hp. Hence, 
EL is renewed as {Ci)?i? 2 }- 

(2-3) Take as a candidate of a head term. According to the algorithm LowCand, 
CL = ELULList = {'Ci)'?iC 2 }- In Example 16, is in Hp. By the algo¬ 

rithm renewJow, LL = {^f}, LList = {^i}, RR = EL\LL = {^i'C 2 } and EL is 
renewed as the empty set. 

(2-4) Take as a candidate of a head term. Then, as Neighbor(LL) = {'Ci)CiC 2 } and 
?2 ^ ^ 1)^2 ^ ?i? 2 , we obtain UU = {^^,^ 1 ^ 2 } and EL = 0 U RR = {^ 1 ^ 2 }- Hence, 
CL = EL ULList = {Ci)?i? 2 }- Since — ^Ci ^2 is in Hp by Example 16, then 
LL = LList = and RR = EL\LL = 0. 

(2-5) Take as a candidate of a head term. Then, as Neighbor(LL) = {CiC 2 ,CiC 2 } 
and ^^^2 ^ CiCl) && -< we obtain D = Moreover, U = 

{Ci;?iC 2 }- UU is renewed as {^i,^i^ 2 }\U = 0 and RR is renewed as t/ U RR = 
As n RR = EL = U RR = e?e 2 ", and CL = 

EL U LList = {^ 1 , ^ 1 ^ 2 , Cl, C 1 C 2 J CiCl}- Since C 1 C 2 can not be a head term in Hp by 
Example 16, the set EL is renewed as ELUjCiCl} = {Ci> CiC 25 CiCl) CiClI- 
(2-6) Take Cl as a candidate of a head term. By the algorithm LowCand, CL = 
ELULList = {C?,CiC 2 ,C|,C?C 2 ,C?C|,CiC|}- In Example 16, C| - yC^Cl + Cl is in 
Hp. By the algorithm renewJow, LL = {C^CiCll- LList = {Ci,C 1 C 2 ,Ci,C 1 C 2 
RR = EL\LL = {CiC 2 ,CiC 2 } and EL is renewed as the empty set. 

(2-7) Take Cl as a candidate of a head term. Then, as Neighbor(LL) = {Cf, CiC 2 ,CiC 2 ! 
C^Cll and CfC 2 /Ci; C^Cl/Cii C^Cl/Ci ^ ht(SList) U MList ULList, hence, by the al¬ 
gorithm LLem, D = {C|}. Moreover, EL = U RR = {C 1 C 2 ,C 1 C 2 .Ci} and CL = 
EL U LList = {C^C?C 2 ,C?C 2 ,C?C 2 ^ClC 2 ^Cl,C|}■ 

3.3. On danger strata 

Here, we present an algorithm for computing bases of algebraic local cohomology 
classes Hp, on danger strata. Basically, we follow the first part of the main algorithm 
ALCohomology for danger strata. However, we can not directly follow it, because the 
termination of the algorithms MonoSafe and BodySafe, is not guaranteed beforehand. 
If (GP) is not a zero-dimensional ideal in (*1) of MonoSafe, then a number of elements 
which do not belong to (G), may not be finite. In such a case, the algorithm does not 
terminate, and this means that (F) is not a zero-dimensional ideal in Ar[[a;]]. In order 
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to resolve this matter, the following algorithm for danger strata is introduced instead of 
the algorithm MonoSafe. The termination of following algorithm is guaranteed by the 
same reason of the algorithm MonoSafe. 


Algorithm 11. (MonoDanger) 

Specification: MonoDanger(A, GP) 

Computing monomial elements of bases of on a danger stratum A. 

Input: (A, GP): a segment of a CCS of F s.t. for all a € A, {aa{F)) is nonzero¬ 
dimensional in A [a:]. 

Output: M : a finite set of triples {A',M,G) where M,G C AT[5], If (G) is zero¬ 
dimensional in K[^, then the set M is MList of Flp on A', otherwise, M = 0. 

BEGIN 

Af ■«— 0; B ^ compute a CCS of Mono(CV(GP)) on A in K[^] 

while P 7 ^ 0 do 

select (A', GP') from B- B ^ S\{(A', GP')}; G ^ ht(GP') 
if dim({G)) = 0 in A'[^] then 

M ■<— compute monomial elements which do not belong to (G) 

M ^ Af U{(A',M,G)} 
else A1 ^ A1U{(A',0,G)} 

end-if 
end-while 
return Af 
END 


The termination of the algorithm BodySafe is also a matter of grave concern on 
danger strata. We have two ideas to resolve this matter. 

The first idea is preparing a natural number u which is an estimated bound of a 
dimension of the vector space Hp.hi many cases, a natural number u can be computed 
from the input P. For instance, if / is a Newton non-degenerate polynomial defining an 
isolated singularity at t he origin G, a bound o f the dimension Hp can be computed by the 
Kouchnirenko formula ( Kouchnirenko . 1976ll . where P = ■ • ■, If a number of 

elements of MList U SList is bigger than and equal to v, then (P) is not zero-dimensional. 
Otherwise, (P) is zero-dimensional in Pr[[a;]]. 

After deciding the number we can compute bases of algebraic local cohomology Hp 
on danger strata as follows. 

We name the same name (BodyDanger) to both Algorithm 12 and 13. By one’s 
strategy, one can select one of them. 


Algorithm 12. (BodyDanger) the first idea 
Specification: BodyDanger(i/, PP, P) 

Computing bases of a vector space F[p on danger strata 

Inpnt: VC: a set of lists ([A, CT, GList, SList, MList, LList, FL, LL, EL, RR, UU]), 
v: a natural number (an estimated bound). 

Output: {S, V): S = SList^, MList^, FLi]} where SList^ UMListi is a basis of Hp 

on Ai, and FL^ is a set of failed candidates for head terms on A^. P is a set of stratum. 
On a stratum of P, P does not satisfy {a € A|/i(a) = • • • = fp{a) = 0} = {Oj. 

BEGIN 
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5^0; 

while VC 7 ^ 0 do 

select £ = [A, CT, GList, SList, MList, LList, FL, LL, EL, RR, UU] from VC 
VC ^ VC\{£} 

iV ■(— a number of elements of SList U MList 

if fc > iV then _ 

(<^1) of BodySafe 

(<^2) of BodySafe 
(<^3) of BodySafe 

(iSi,X>i) ^ BodyDanger(i/,7^,F); V-^V\j{Vi} 

else I? <— X> U {A} 

end-if 
while-end 
return (5, V) 

END 

The second ideal is the following. Let m be the maximal ideal at the origin O (i.e., 
m = {xi,..., Xn)) and be a sufficient big positive integer. Then, F^(u) = (F) + is 
an ideal supported at the origin O where = (x^', X 2 , ■ ■ ■, xj^). Therefore, bases of the 
vector space Hp can be computed by the algorithm BodySafe. li Hp = Hp 
on a stratum A, it is obvious that (F) = on A. That is, if there exists G N such 
that Hp = Hp on A, then (F) is zero-dimensional on A. li Hp ^ Hp 
on A, there exist some local cohomology classes in a basis of Hp such that the local 
cohomology classes do not belong to Hp . By analyzing such local cohomology classes, 
we can easily guess and prove that Hp has infinite many (systematic) elements which 
are linearly independent, on A. That is, in this case, (F) is not zero-dimensional on A. 


Algorithm 13. (BodyDanger) the second idea 
Specification: BodyDanger VC, {/i,..., fp}) 

Computing bases of a vector space Hp on danger strata 

Input: VC: a set of lists ([A, CT, GList, SList, MList, LList, FL, LL, EL, RR, UU]), 
p: a natural number (a sufficient big number). 

Output: {S,V): S = lJj{[Aj, SList^, MList^, ELj]} where SList^ UMListi is a basis of Hp 
on Ai, and EL^ is a set of failed candidates for head terms on A^. F is a set of stratum. 
On a stratum of V, F does not satisfy {a G A|/i(a) = • • • = fp{a) = 0} = {O}. 

BEGIN 


H ^ F^i„) ^ {/i,...,/p,x5',x^,...,x(;};F,^(.+i) ^ {/i,...,/p, xJ'+S x^+\ . 

HiSafeBody(F£, F,p(,.)) ; H2 ^SafeBody(F£, F,^(„+i)) 

S i — Hi n H2i Fl i — Hl\Sl 

while Fi 7^ 0 do 

select £ = [A^, SList^ MList\ FL^j from Vi] Fi Vi\{£'\ ; V ^ Vd {A^} 

end-while 
return {S, V) 

END 




We illustrate the second idea with the following example. 
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Example 21. Let us consider Example 4, again. The term order is the total degree 
lexicographic term order such that ^ ^ 2 - If the parameter t belongs to V(t —2) or V(< + 
2), the ideal (F) is not zero-dimensional in K[x]. Set i/ = 4 and ^^( 4 ) = {-§^, 

and E,„( 5 ) = We apply the algorithm BodySafe for {[V(t — 2),CT, 

Cl el, e?}}, SList, MList, LList, FL, LL, EL, RR, UU]} with , 

where CT, SList, MList, LList, FL,LL, EL,RR,UU are empty sets. Then, a set 
Gi = {1, 6 , Cl, el, 6 ^ 2 , el, - Cl,e? - 6 ^ 2 , ^16 -6^1} fs a basis of the vector space 

^ ^2 = {i,e 2 ,ei,el,eie 2 ,ei,eie 2 - e 2 ,ei - eiel,eie 2 - eieljei - elel+ell 

is a basis of the vector space Flp ( 5 , and Gi ^ G^. 

Our implementation can compute bases oi FIp , . and oi FIp , ,, too. One can find 
some systematic elements based on a certain rule from these bases, and guess that the 
following four sets 

{E-io(-irei”'*elifc = 2n + 4,n e N}, {E-io(-irei”'“e|*^'l^ = 2n + 4,n G N}, 

{E'io(-l)*ei“'*+'elifc = 2n + 4,n G N}, {E■io(-l)*ei"-'*+'el*+'1^ = 2n + 4,n G N}, 
are included in a basis oi Hp on V(t — 2). This can be easily proved. Therefore, (F) is 
not zero-dimensional on V(t — 2). One can also easily verify the non zero-dimensionality 
of {F) on V(t -I- 2). 

In Fig. 6 , we represent an monomial element of Hp as • and an elements of head terms 
of the systematic elements as A. 
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• • 


( 0 , 0 ) 


A A A A . 

. A A A A 
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Example 22. Let fi = x\ + X 2 + sx|x 3 -I- tx 2 X §,/2 = X 2 -I- Xg. It is described in 
( Aleksandrovl 1 1 98^ that /g = /2 = 0 defines a quasi-homogeneous complete intersection 
isolated singularity provided that the parameters s, t do not belong to V((s-|-t)^-|-(s-|-l)^)) 
and the Milnor number is equal to 16. 

Let /g = 3 x 2 X§ -|- 2 SX 2 X 3 -I- &x| — sx| — 21 x 2 X 3 , f 4 = xix §,/5 = xiXt and set 


F = I fi ■ 0. fg. fi. h\. Since /i ,/2 are quasi-homogeneous, a result of (IGreuel. G.-M, 


1975 ) on Milnor number and the Grothendieck local duality theorem (jGrothendieck 


19671) imply that 77j,’ is a vector space of dimension 16 provided /i = /2 = 0 has 


an isolated singularity at the origin. However, the algorithm ZeroDimension outputs 
V(s — t — 1)\ V(t^ -I- 2t^ -I- 2t -I- 1, s — t — 1) as a danger stratum and BodyDanger (our 
implementation) judges {a G X|/i(a) = /2(a) = • • • = /5(a) = 0 } 7 ^ {O} on the stratum. 
One can check the fact V(s — t — 1)\V(t^ -I- 2t^ -|- 2t-|- 1, s — t — 1) ^ V((s -|-1)^ -I- (s -I- 1)^). 
For instance, take (s, t) = (i, —i) G V(s — t — 1)\ V(t^ -|- 2t^ -|- 27 -|- 1, s — t — 1), then 
ih-h) i V((s-k 7)3-h (s-h 1)3) and {a G X|cr( 1 _ i, (/i)(a) = ••• = cr(i _ i)(/ 5 )(a) = 
0} 7^ {O}. 

The algorithm BodyDanger works powerfully to find strata on which {a G X|/i(o) = 
• • • = /p(a) = 0} 7^ {O}. 
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We conclude this section by briefly discussing the effectiveness of the proposed method. 
In order to detect unnecessary strata on which (F) is not zero-dimensional in first, 

the algorithm ZeroDimension decomposes the parameter space into safe strata and 
danger strata. If there exist danger strata, second, the algorithm MonoDanger detects 
unnecessary strata. After that if there still exist undeterminable strata, finally, the algo¬ 
rithm BodyDanger detects unnecessary strata by computing local cohomology classes 
ot Hp on the strata. The final step is a practical method for detecting unnecessary strata. 
Note that if we compute a parametric local cohomology system without the algorithm 
ZeroDimension, then the algorithm BodyDanger (the general case is the second idea) 
has to be always performed because all strata of are regarded as danger. As we de¬ 
scribed above, BodyDanger actually computes local cohomology classes several times. 
Thus, in this case, the computational complexity increases. To avoid an increase in com¬ 
putation cost, we have innovated the algorithm ZeroDimension. As we described in 
Example 22, the algorithm ZeroDimension powerfully helps for checking unnecessary 
strata, and makes the computational method of a parametric local cohomology system, 
more effective in computational speed and complexity. 


4. Parametric standard bases 

Here, we introduce an algorithm for computing parametric standard bases of zero¬ 
dimensional ideals by using bases of algebraic local cohomology classes. 

Definition 23 (inverse orders). Let ^ be a local or global term order. Then, the inverse 
order of -< is defined by x°‘ -< x^ a:“. 

If -< is a global term order (1 is the minimal term), then is the local term order 
(I is the maximal term). Conversely, if -< is a local term order, then is the global 
term order. 


4 . 1 . Parametric standard bases 

Definition 24. Let F be a subset of (Ar[t])[[x]], a stratum in A'™, Si a subset 
of Ar[t]Aj[x]] and ^ a local term order where 1 < i < /. A finite set S = {(Ai,5'i), 
..., (A;, Si)} of pairs is called a parametric standard basis on Ai U • • • U A; for {F) 
w.r.t. -< if aa(Si) is a standard basis of the ideal (aa(F)) in .A[[x]] w.r.t. ^ for each 
i = 1,..., I and d € Ai. 


Let F = {/i,...,/p} be a set of polynomials in (A'[t])[x] such that generically 
{a G A| /i(a) = • • • = fp(a) = 0} = {O} where A is a neighborhood of the origin O 
of AT". Then, by utilizing the information of bases oi FI p, one can obtain parametric 
standard bases of (F) in Ar(t)[[x]]. 

Let us recall that there is a natural pairing, denote by reS{o}( i )i between the quotient 
space Ar[[x]]/(F) and the vector space iJp where (F) C Ar[[x]] is a zero-dimensional ideal. 


reS{o}( , ) : K[[x\]/{P) x Hp 


K. 


Since the pairing is n on-degenerate according to the Grothendieck local duality theorem 
(IGrothendieckl . Il957l) . we have the following result. 
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Lemma 25 ( Taiima and Nakamural (120091) '). Let P = {gi, ■ ■ ■, gq} be a set of poly¬ 
nomials in K[x\ such that {a G X\ gi{a) = ••• = gq{a) = 0} = {O}. Then, a given 
formal power series h G [[a;]] is in the ideal (P) if and only if for all G h satisfies 
reS{o}(/i, ip) = 0 where dt is a basis of the vector space Hp. 

One can extend this fact into the parametric cases. The next theorem gives us the 
relation between bases of Hp and parametric standard bases of (F). 


Notation 26. Let SList be a set of polynomials in and LList be a set of lower 

terms of all elements of SList. Suppose that there is no monomial in SList, and SList 
has an element whose form is -I- C(a,k)C^ where G K{t). Then, the transfer 

^^(SList,LList) is dehned by the following: 

f-S' 5(SList,LList)(?^) = ^ in K{t)[x\ if G LList, 

\ Ght(SList) 

I ^^(SList,LList) in K{t)[x\ if ^ LList. 

Let G be a set of terms in K[^]. Then, the set ^^(SList,LList) (G) is also defined by 

5i?(SList,LList) (G) = {*S'i?(SList,LList) & G}. 

Theorem 27. Let ^ be a global total degree lexicographic term order (Definition 1). 
Let (5,2?) be an output of ALCohomology(F) and a list [A, SList, MList, LList, FL] 
is in S. Then, for all o G A, cra(523(SList,LList)(FL)) is the reduced standard basis of 
{aa{F)) w.r.t. (the local total degree lexicographic term order), in .A[[a;]]. Namely, 
{(A, 5i2(SList,LList)(FL))} is a parametric standard basis on A for F. (The notation a is 
from section 3.1.) 


Proof. Since the algorithm BodySafe decides linear combination elements of a basis 
of Hp from bottom to up w.r.t. ^ and (F) is zero-dimensional on A, the set CV~^(FL) 
(failed candidates of head terms) becomes a set of head terms of the standard basis w.r.t. 

on A. By Lemma 25 (an d Theorem 7, Proposition 8 and Theorem 9 of the paper 
( Taiima and Nakamural[20091) ). for all a G A, it is obvious that if G FL is not in LList, 
then the monomial itself is in the ideal (cra(F)) in .A[[x]], and if G FL is in LList, 

then ffa \ x^ — ^ is also in {aa{F)) and is not in ht(SList) w.r.t. 

\ ^''Ght(SList) J 

Hence, for all a G A, erg (5i?(SList. LList) (FL)) is the reduced standard basis of {aa{F)) 
w.r.t. on A. □ 


This theorem leads us to construct the following algorithm for computing parametric 
standard bases. 


Algorithm 14. (StandardBasesl) 

Specification: StandardBasesl (F) 

Computing a parametric standard basis for a zero-dimensional ideal (F). 

Input: F G (A[t])[x], a global total degree lexicographic term order. 

Output: (5, A 2 ): 5 is a set of pairs (A, E) such that for all a G A, <Ja{E) is the reduced 
standard basis of {aa{E)) w.r.t. A 2 is described in the algorithm ALCohomology. 
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BEGIN 

5^0; (^ 1 ,^ 2 ) ALCohomology(E) 
while Ai ^ % do 

select B = [A', SList, MList, LList, FL] from Ai] Ai -ir- Ai\{B} 

S^SU {(A',5B(SList,LList)(FL))} 

end-while 
ret urn (5, A 2 ) 

END _ 

This algorithm gives a nice decomposition of the parameter space depending on struc¬ 
ture of bases oi Hp thanks to the algorithm ALCohomology. This is the big advantage 
of this algorithm. 

Example 28. Let / = be a polynomial with a parameter t in (C[t])[a:i, ^ 2 ] 

and ^ be the global total degree lexicographic term order such that -< ^ 2 - Set F = 
■§^}- The output of ALCohomology (F), is already given in Example 16. 

(T) On V(t), SList = 0,MList = {1, 6 , 6 ,, 66 ,Cl,LList = 0 and 

FL = { 6 , 6 n- 

(ii) On C" \ V(t (6 - 4)), SList = { 6 ?| - £.1 - UU 2 AI - HKI + ?t}, MList = 

{i, 6 , 6 ,e?, 66 , e|}, LList = and FL = { 6 ,Cl}-_ 

In case (i), 5i?(gList LList)(FL) = {xfjxl} is the reduced standard basis of (F) w.r.t. ^ 

In case (ii), each elements of FL is transformed as follows: 

6 S' xl + ^xix"^, 66 —S' xlx2 + jxl, 6C| —s> xixl, £l —)• x^. 

Therefore, {xf + ^xiX 2 ,xlx 2 + fxl,a;fx 2 ,a;f} is the reduced standard basis of (F) 
w.r.t on C \ V(t (6 — 4)). 

Let us remark that if t = ±2, then (F) is not zero-dimensional in F[[a:]]. 

All algorithms of this paper have been implemented in the computer algebra system 
Risa/Asir by the authors. In the following example, we give an output of our implemen¬ 
tation. 

Example 29. Let F = {3sxl + 2xiX2+tx2,2xiX2+5x2+3txiX2} be a set of polynomials 
with parameters s,t in (C[s,t])[a;i, 0 : 2 ], and ^ be the global total degree lexicographic 
term order such that xi -< X 2 - Generically, F has only the point O in A where A is a 
neighborhood of the origin O of C^. The variables 6,6 s-re corresponding to variables 
a;i,a; 2 . Our implementation outputs bases of the vector space Hp and standard bases of 
(F) on w.r.t as follows. 

- If the parameters belong to V(s), then (F) is not zero-dimensional. 

- If the parameters belong to \ V(st(—15s-f 2t)), then a set {1, 6 , 6 , ■fl, 66 , 

66 , + § 566^2 + is a basis of Hp. Hence, 

{s^tx1 + {—^s+^t)x 2 + ^st‘^X 2 , s‘^txiX 2 + {^s^ — ^st)x 2 , X 2 } Is & parametric standard 
basis w.r.t 

- If the parameters belong to V(—15s + 2<)\ V(s, t), then a set {l, 6 , 6 ,^ 2 , 66 , s 6 ~ 
6 ^ 1’*^2 “ 6 ^ 16 } is a basis of Flp. Hence, {sxf + 62 ^ 2 ,^ 1 ^ 2 ,^ 2 } is a parametric 
standard basis w.r.t 

- If the parameters belong to V(t)\V(s,t), then a set {1, 6 , 6 , ^ 2 , 66 , , ^6^1 ~ |Ci, 

6^2 + s 6^2 ~ i?i 6 } is a basis of Hp. Hence, {sxf + ^xiX 2 , —^xix^ -f sxj} is a 
parametric standard basis w.r.t 
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Note that in case Y{t)\Y{s,t), the dimension of the vector space Hp is 8, but in cases 

\ V(st(—15s + 2t)) and V(—15s + 2t)\V(s, t), the dimension of the vector spaces Hp 
are 7. Our implementation tells us this deference. 

4-2. Other local term orders 

The algorithm ALCohomology has been constructed based on a global total de¬ 
gree lexicographic term order -<i. That’s why reduced standard bases w.r.t. are 

directly obtained by outputs of the algorithm ALCohomology. Here, we describe how 
to compute standard base w.r.t. other local term orders. 

Let (*S,T’) be an output of ALCohomology(J^) and ^ be a global term order in 
K[^]. Suppose that [A, SList,MList,LList,FL] G S, SList = {-ipi, ■ ■ ■ ,ipp} C K[^] and 
the list is lined up all elements of Mono(SList) in order of -< where 

^ ^ ... ^ Moreover, let M the coefficient matrix of SList w.r.t. the vector 

(i-e-, = M where , V’p) is 

the transposed matrix of {ipi,... jipp). Then, it is possible to compute the row reduced 
echelon matrix of M on A, like a method for solving the system of parametric linear 
equations. Let {(Ai, M{), (A 2 , M 2 ),..., (A/, M/)} be a set of pairs such that for each 
1 < i < I, Ml is the row reduced echelon matrix of M on A^ and A = Ai U • ■ • U A;. Then, 
we have the following theorem. 

Theorem 30. Using the same notation as in above discussion, let *((^ 1 ,(^ 2 , • ■ • ,‘Pp) = 
M' where \ < i < 1. Suppose that SL = {(^ 1 ,(^ 2 ,-•• , ‘Pp}i L = 

Mono(S'L)\ ht(S'L), T = ht(S'L)UMList and the reduced Grobner basis of (Neighbor(r) 
\T) is FL^ w.r.t. -<. Then, for all a G Ay aa{SB(^sL,L){FL^ j) is the reduced standard 
basis of {aa{F)) w.r.t. in ^"[[ 0 :]]. (The transfer SB(^sl,l) is from Notation 26.) 

Proof. As SList U MList is a basis oi Hp on A, it is obvious that SL U MList is a basis 
of Hp on Ay too. L is a set of lower terms of SL w.r.t Since M/ is the row reduced 
echelon matrix of M w.r.t. the vector *(C“U)•■•)?“’') on A, the set CV~^(FL^) 
becomes a set of head terms of the standard basis w.r.t. on Ay By this observation 
and Theorem 27, this theorem holds. □ 

This theorem leads us to construct the following algorithm for computing parametric 
standard bases w.r.t. any local term order. 


Algorithm 15. (StandardBases2) 

Specification: StandardBases2(F, -<) 

Computing a parametric standard basis for {F) w.r.t. 

Input: F C (Ar[t])[a;], a local term order. 

Output: [S, A 2 ): 5 is a set of pairs (A, E) such that for all a G A, <Ja{E) is the reduced 
standard basis of {aa{E)) w.r.t. A 2 is described in the algorithm ALCohomology. 
BEGIN 

5 -G- 0; (Ml, A 2 ) ALCohomology(F) 
while Ml ^ 0 do 

select B = [A', SList, MList, LList, FA] from Ai; Ai Ai\{S} 

V <— Line up all elements of Mono(SList) w.r.t. 
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M <— Make the coefficient matrix of SList w.r.t. v 

AA4 Compute the row reduced echelon matrix of M on A' 

while AA4 ^ 0 do 

select (A", M') from AM] AM ^ M')} 

(where M' is the row reduced echelon matrix of M on A".) 

\aiM 2, ■■■ ,(Pp) <- M' *V] SL ^ {(/?!, (/? 2 , • • • ,Ap} 

L -s— Mono(S'L)\ ht(S'L); T -s— ht(S'L) U MList 
FL^ ^ the reduced Grobner basis of (Neighbor(r)\T) 
5^5U{(A",5B(Si,i)(FL^))} 
end-while 
end-while 
ret urn (iS, A 2 ) 

END 

Example 31. Let / = xf+tx^x'^+x^ be a polynomial with a parameter t in (C[f])[a;i, X 2 ]. 
Set F = The output of ALCohomology(F), is already given in Example 16. 

If a parameter t belongs to -4)), then {1,Ci,6 ,,Ci? 2 , ^1,6^2 “ “ 

^?i? 2 ,C 2 “ i?i ?2 + Cl} is a basis of Hp- Let ^ be the local lexicographic term order 
such that xi -< X 2 and SL = “ yCiCl +Ci}- 

We compute a parametric standard basis of (F) w.r.t. ^ on C\V(t(t^ — 4)). First, 
we line up all elements of Mono(SList) w.r.t. (the global lexicographic term order), 
then we get the vector v = *(Ci > Ci > C1C2 1 C1C2,CiCl^Cl i Cl)- The coefficient matrix of 
SList w.r.t. V is M, and the row reduced echelon matrix of M is M': 

Cl C? C?Cl CIC2 CiCl cl Cl 

/o -l/2t 0 0 1 0 o\ /lo-i/to 0 1 0 \ 

M = 0 0 0 -l/ 2 t 0 0 1 , M' = 0 1 0 0 - 2 t 0 0 . 

\1 0 -1/t 0 0 10/ \0 0 0 1 0 0 -2t/ 

Then, M' ‘u = ‘(C| - iClC2^C? - 2 ^ClC 2 ^C?C 2 - 2 te|). Hence, SL = - iClC2^Cl - 

2^ClC2^C?C2-2^Cl}, ht(FL) = {Cl,C?,C?C2}, L = Mono{SL)\U{SL) = ClC2^Cl,Cl} 
and T = ht(S'L) UMList. Since the set FL^ = {Ci? C1C2) CiCIj ClI is the reduced Grobner 
basis of (Neighbor(r)\r), the parametric standard basis of (F) w.r.t. ^ on C\V(t(t^ — 
4)) is 

SB(^SL,l){FL^) = {x\,x\x2,xixl + 2txl,xl + 2txlx2}- 
In Fig. 7, we represent an element of MList as •, an element of ht(S'A) as A. As the set 
FL^ plays a key role to construct standard bases, we give the elements of FL^ in the 
figure. 



Fig. 7 
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5. Conclusions 


We have described algorithms for computing parametric local cohomology systems, 
and given a new algorithm for computing parametric standard bases as an application. 
The algorithm for computing parametric standard bases, has the following advantages. 

- The algorithm always outputs a “reduced” standard basis. The computer algebra 
system Singular has a command that outputs a (non-parametric) standard basis. 
Singular does not have this property. 

- The substantial computation consists of only linear algebra computatio n. 

- We do not need Mora’s reduction (tangent cone algorithm (jMora . 1982h l for com¬ 
puting standard bases. 

- The algorithm outputs a nice decomposition of the parameter space depending on 
the structure of standard bases w.r.t. a local total degree lexicographic term order. 

All algorithms of this paper, have been implemented in the computer algebra system 
Risa/Asir. Actually, there does not exist any implementation for computing “paramet¬ 
ric” standard bases, except for our implementation. Only our implementation exists for 
them. Our implementation is useful for studying and analyzing singularities. 
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